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O ' Abstract 



The ground-state properties of the S = 1 Haldane-Shastry model are studied 
using a modified Lanczos algorithm and diagonalizing exactly small chains. We find 
evidence that, as for the antiferromagnetic Heisenberg model, the spectrum shows a 



O . 1 

O ■ gap, in contrast to the S = ^ case. The correlation functions < S^{0)S^{'m) > decay 



exponentially for large m. We find that the correlation functions for the Haldane- 
Shastry model decay faster than for the Heisenberg model. We estimate the infinite 
system limit for the ground-state energy, value of the gap and correlation functions. 
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The S = J Haldane-Shastry model [1,2] has attracted considerable attention since it 
is an integrable model that belongs to a class of systems with long-range interactions. 
The ground-state energy and correlation functions have been obtained [1-3] together with 
the thermodynamics [4]. The ground-state wave function is a spin singlet of the Jastrow- 
Gutzwiller form. The excitations are spin-| spinous [4] that form a gas of a semionic 
nature [4,5]. The asymptotic correlations decay algebraically with exponent rj = 1 without 
logarithmic corrections, in contrast to the Heisenberg case. This indicates the absence of 
spin exchange between the spinous rendering the model solvable in greater detail than the 
short-range Heisenberg counterpart, solvable by the traditional Bethe ansatz method. 

On the other hand, the similarities to the Heisenberg model lead to the interest in 
studying the model for higher values of the spin. In particular, the case of spin S* = 1 
(and in general the integer spin cases) are of interest since it has been established on 
general grounds and verified explicitly for short-range interactions, that these values are 
qualitatively different from the half-integer ones due to the Haldane gap [6]. 

The model Hamiltonian is a periodic version of -^ exchange given by 

i<j sm [(j){t-j)] 

where 4> = ^ with A^ the number of lattice spins with lattice constant 1. We consider 
only the antiferromagnetic case J = 1. 

In this paper we consider the S* = 1 Haldane-Shastry model (HS) and study the 
ground-state properties. We use a modified Lanczos method [7,8] to obtain the ground- 
state and first excited state energies and the correlation functions for small systems of size 



up to 16 spins. We compare our results to those previously obtained for the Heisenberg 
case stressing the similarities and differences. We estimate the large system size limit of 
these quantities using the extrapolation method of Vanden Broeck and Schwartz (VBS) 
[9-11]. 

The exact energy spectrum has been obtained for the S = ^ case [1-4]. The groundstate 
is a spin singlet with energy [2] 

E n^ / 5 \ 

ivj = -24 (' + ivj) '2' 

with a -^ correction like for the Heisenberg model [11]. The first excited state for A^ even 
and for the Heisenberg model is a S't = 1 triplet where St is the total spin of a A^-spin 
state. In the case of the HS model the first excited state is however a S't = 1 quartet. 
(For A^ odd the groundstate is a quartet for both models). The full spectrum has been 
obtained [4] showing a supermultiplet structure such that the energies in units of | f -^j 
are integers [1]. A set of states with A''^^ spin-| spinon excitations with St = S^ = -^ 
(with M = {N — Nsp)/2 reversed spins) generates the full energy spectrum [4]. 

The Hamiltonian eq. (1) can be exactly diagonalized if we consider small chains [1]. 
For spin-S* we have a basis of (25* +1)^ states which grows very fast as S and A^ grow. 
Taking A^ < 8 we have obtained the full spectrum for S* = 1. For A^ even, the groundstate 
is again a singlet and the first excited state is a 5*^ = 1 triplet, as for the Heisenberg case. 
For A^ odd the groundstate is a singlet in both models. 

Since the size of the matrices to diagonalize increases exponentially with the size of 
the chain, we used the modified Lanczos method [7] to obtain the groundstate and the 



first excited state energies, together with the groundstate correlation functions. 

The size of the vectors can be considerably reduced using the symmetries of the prob- 
lem. The Hamiltonian commutes with the total spin operators {St) and S^-, with the 
translation operator T, the spin flip operator R and the reflection operator L {i ^ 
N + 1 — i,i = 1, . . . , A^). The ground state has total 3^ = and one of the (degen- 
erate) first excited states also has 3^ = 0- One can then immediately reduce the states 
under consideration to this subspace only. 

The operators T and L commute with R, i.e. one has [T, R] = [L, R] = 0, but they 
do not commute with each other, since LT = T^^L. The ground state eigenvalues are 
r = I = 1 and t = 1, corresponding to an wavevector A; = 0. The first excited state has 
r = I = —1 and t = — 1, corresponding to A; = vr (for an even number of lattice points). 

Although in general the operator L mixes k and —k states, the ground state and the 
first excited state under consideration are simultaneous eigenvectors of the operators T, R 
and L. 

The action of the local operators 3^ and S'f is simply given in the direct product 
basis \ml > ■ ■ ■ |m^ >. In general, these states are not eigenvectors of these additional 
symmetries. One forms then the classes of states which are closed under them. One 
starts with a state \a >, applies the translation operator T rit times until one finds 
T"'|a >= I a >, where rit is necessarily a divisor of A^. One proceeds similarly with the 
operators R and L for which rij. = 1,2 and rii = 1,2. The multiplicity of this class 
is ntUrni. The state |a > is then the representative of this class. The other classes 
are formed proceeding in the same manner starting with other states, not already used. 



until the state space is exhausted. An eigenvector of the operator T (with eigenvalue t 
obeying to t"* = 1) is obtained applying to its representatives the symmetrization operator 
n[ = (l+t~^T+(t~iT)2+...+ (t-iT)"*-i). For the states in this class we have f]; = ^^t, 

where the general operator l^^ = (1 + t'^T + (t^^T)^ H h (t-^T)^'^) and St{a) = N/rit 

is the symmetry factor. The same applies to the operators R and L. As stated above, 
when t = ±1 it is possible to construct a simultaneous eigenvector of the operators T,R 
and L. In the cases in which the operator R or L does not introduce new states (n^ = 1 
or fir = 1) one should verify that the eigenvalues induced by the first operator (s) are 
compatible with those of the other operator (s). This allows to reduce even further the 
number of classes. Any eigenvector with these eigenvalues is completely defined by its 
projections on the representatives of the classes. 

The method goes as follows [7] . To obtain an approximate ground state wave function 
we choose a trial wave function ipQ that can not be orthogonal to the true ground-state. 
We define a state ipi as [7] 

HiJo- <H>iJo 
{< H^ > - < H >^) ' 

where < '?/'i|'?/'i >= 1, < ipi\ipo >= and < if" >=< ipolH^lipo >■ Defining a matrix 
of the Hamiltonian in the basis ipo^'^i we can diagonalize it [7] obtaining a next order 
approximation for the energy Ci, and ground state wave- function ipQ, with 



e^ = < H > +ba (4) 



^0 = .. . ,.1/0 (5) 



where 

b = [<H''> - <H >')''' (6) 

« = f-{f + lf" (7) 

and 

_ <i/3 > -3 < H >< H^ > +2 < H >^ 



2 (< /J2 > _ < /J >2)3/2 

Taking ipo as the new ipo we can iterate the method to obtain a better estimate for the 
energy and ground state wave function. 

For the first excited state we proceed in a similar manner, starting from a trivial wave 
function orthogonal to the true ground state (as guaranteed by some symmetry argument), 
but not orthogonal to the true first excited state. For small chains, where the use of the 
symmetries is not strictly necessary, one can look for the excited state with S^ = 1, which 
is necessarily orthogonal to the ground state. However, for larger chains one has to take 
into account the different symmetries, and with the additional eigenvalues one can still 
look for the S^ = first excited state, as explained above. 

Since the complete Hamiltonian commutes with the operators T, R and L only transi- 
tions to the classes constructed above are allowed, even if separate terms allow them. The 
action of the Hamiltonian on a state is obtained writing, within each class, HQ'\a >= 
-4^QH\a >, where Q is the product of the three symmetry operators for T, R, and L. 
If H\a >= J2b o^L^' R^''TP^\b > , where pi^r,t are integers and \b > is the representative of 
a class, one finally finds < b\Hi^'\a >= Z);, f74Q;/^'r'''"tP*. One also has to take into ac- 
count the multiplicity of the classes when normalizing states and making inner products 
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of states. 

The accuracy of the method was tested against the exact diagonahzation of the full 
spectrum for sizes up to iV = 8. We considered system sizes up to A^ < 16 (with N even). 

This procedure gives the results for the several finite-size systems. We are however 
interested in the infinite system limit and standard extrapolation methods [9] have to be 
used like the VBS method [9,10]. In this method we want to estimate the limit of a finite 
sequence P„ (n = 1, ..., A^). Defining 

p(m+l) _ p{m.) _| (q\ 



1 



-* n+1 ~ -< « 

where Q^"^-* = 0, P^°^ = P„ we obtain an estimate of the sequence iterating (like the 
recursive rule due to Wynn). If am = this is the Aitken-Shanks transformation which 
is adequate for exponential behavior. To generate the Pade-Shanks transformation we 
select am = 1- A power law behavior is well fitted choosing the Hamer and Barber's 
transformation am = — [1 — (— 1)'"]/2. We get an estimate of the asymptotic value of the 
sequence P„ [9,10] selecting am appropriately. 

i) Ground-State Energy 

In Table I we show the values for the groundstate energy of the S = 1 Haldane-Shastry 
chain with sizes A^ = 4 to 16 (A^ even). A linear fit of the groundstate energy per spin as 
a function of -^ yields (1.2568 ± 0.0016). Using the more accurate VBS method yields 
the results shown in Table II. These results indicate that the behavior is not purely -^ in 
contrast to the S = ■^ case. Other terms have to be included, particularly for small values 
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of A^. Since the spectrum has a gap (see below) we expect that the exponential behavior 
is more adequate [12,13]. It has been found for the Heisenberg model [13] that am = 1 
yields the best results. This has been found looking at the decay length at size N and 
requiring that 

e(n,m)= ) (11) 

ni I p(m-l)_p(m-l) I 

be such that 

^{n,m) <^{n + l,m-l). (12) 

Only the Pade-Shanks transformation yielded consistent results [13]. In our case, however, 
all three cases yielded consistent tables and therefore we include the three choices. 

ii) Gap 

We consider now the gap (the difference in energy between the first excited state and 
the groundstate). In Table III we give the values of the gap as a function of A^ for the 
Haldane-Shastry model. In Fig. 1 we show the gap and the gap per spin as a function of 
A^ for the two models. The results seem to indicate that the gap is finite for both models 
in the A^ — > oo limit in agreement with Haldane's conjecture [6] and with previous results 
for the Heisenberg case. The VBS extrapolated values are shown in Table IV. In the case 
of the gap the Pade-Shanks yields once again consistent results but the other two do not 
obey eq. (12) for the last iteration of eqs. (9,10). We recall that, for comparison, the 
value of the gap for the Heisenberg model is 0.41050 [13,14]. The extrapolated value of 
the gap per spin is close to but not zero due to the small number of points considered. 



iii) Correlation Functions 

The groundstate correlation functions are defined by 

o N ^ czcz ^ 

In Table V we show them for the same set of system sizes and in Table VI the VBS 
extrapolated values. 

In Fig. 2 we compare the correlation functions with those for the Heisenberg model. 
For a finite system the energy scale of the two models is not the same. However, in the 
infinite system limit the nearest-neighbor interaction is the same. We feel therefore that 
it is worthwhile to compare the behavior for the two models. In Fig. 2a |Ci| is displayed 
as a function of A^; in Fig. 2b we show C2 also as a function of A^ and in Fig. 2c we show 
Cm for A^ = 16 as a function of m. The general trend is that, in spite of the long-range 
nature of the interaction, the correlation functions for the Haldane-Shastry model decay 
faster (in the sense that the numerical values are smaller) than those for the Heisenberg 
model. A possible interpretation is that, similarly to the S = ■^ case, there may be an 
absence of interactions between the spinous. 

It has been argued that the behavior of the correlation functions Cn/2 when N -^ 00 
reflects their behavior with distance in the infinite system [15,8]. If there is a gap in the 
spectrum it is expected that 

lim \Cm\ ~ Be^f (14) 

m— »oo 

instead of the power-law behavior observed in the S = ■^ case (in both models). In Fig. 3 



we plot In |CAr/2| vs. A^ for both models. A linear fit yields for the Haldane-Shastry model 

In \Cn/2\ = (-0.02 ± 0.02) - (0.162 ± 0.0019)Ar (15) 

which corresponds to a correlation length of the order of .^ = 3.1. A similar analysis for 
the Heisenberg model gives C,h of the order of 4.9. The quality of the fits is similar for both 
models. However, the correlation length for the Heisenberg model has been estimated to 
be C,H = 6.03 [14] obtained from a fit to the large distance limit of the Bessel function 
Ko(r/^), which asymptotically behaves as ~ {^/rY'^e~^'^, and, therefore, a plot like in 
Fig. 3 underestimates the correlation length for the Heisenberg model. 

In Fig. 4a we show In |C(m)| as a function of m for A^ = 16. At large m we expect 
an exponential decay of the form of eq. (14) with possible corrections. It is clear that the 
exponential behavior has not been fully reached. The exponential behavior is probably, 
as for the Heisenberg model, only to leading order. We have tried several functional forms 
like those for the Heisenberg model. In Fig. 4b we plot F{m) = ln(|C(Tn)|m^/^) as a 
function of m for A^ = 16. It is clear that for these system sizes the fit is poor even for 
the Heisenberg model. One requires larger systems to reach the limiting behavior. Note, 
however, that the functional form for the Haldane-Shastry has not been determined and 
therefore the numerical results are not conclusive to extrapolate to the large A^ limit. On 
the other hand, the linear fit of Fig. 3 is quite good for both models. We take the value 
obtained for ^ for the Haldane-Shastry model to be the estimate obtained by this fit. 
Larger system sizes and the appropriate functional form are therefore required to more 
accurately determine the correlation length for the Haldane-Shastry model. 
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In summary, in this work we have studied the groundstate of the S* = 1 Haldane- 
Shastry model, which has been recently presented [1,2] and solved exactly for S = ^. The 
importance of the model lies on the fact that it belongs to a new class [5] of integrable 
systems. We have confirmed that, according to arguments by Haldane, the model shows 
a finite gap in the energy spectrum. We based our result on the extrapolated value of the 
calculated gap for various small chains and on the exponential behavior of the correlation 
function Civ/2 with A^. We have also obtained, by comparison to the Heisenberg case, that 
the correlation functions for the Haldane- Shastry model decay faster as a function of both 
m and N. The same happens in the S = 1/2 case where the spectrum is gapless. This 
is due to the oscillatory nature of the (positive) interaction in this model. The estimated 
correlation length is smaller for the Haldane-Shastry model and the value of the gap is 
larger than for the Heisenberg model. These exact numerical results further extend the 
similarities and differences between the two models for a case {S = 1) which is not exactly 
solvable, in contrast to the 5* = | case. 
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TABLE I- Groundstate energy and groundstate energy per spin as a function of A^. 

N -En -En/N 

4 6.168503 1.542126 

6 8.270682 1.378447 

8 10.59553 1.324441 

10 13.00953 1.300953 

12 15.46690 1.288908 

14 17.94784 1.281988 

16 20.44263 1.277664 



Table II- Extrapolated values for —En/N for the Hamer-Barber (H-B), Pade-Shanks 
I^P-S) and Aitken-Shanks (A-S) transformations eqs. (9,10). 

H-B P-S A-S 
1.263147 1.267894 1.265328 
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TABLE III- Gap between the groundstate and the first excited state as a function of 

N. 

N Gap Gap/N 

4 1.2337 0.30843 

6 0.87226 0.14538 

8 0.71962 0.08995 

10 0.64551 0.06455 

12 0.60678 0.05057 

14 0.58550 0.04182 

16 0.57333 0.03583 



Table IV- Extrapolated values for the gap for the Hamer-Barber (H-B), Pade-Shanks 
(P-S) and Aitken-Shanks (A-S) transformations eqs. (9,10). 

H-B P-S A-S 
0.55045 0.55439 0.55285 
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TABLE V- Correlation functions Cm for the set of values A^ = 4 to 16 with A^ even. 

N 4 6 8 10 12 14 16 
m 

1 -0.75000 -0.71660 -0.70630 -0.70210 -0.70017 -0.69923 -0.69876 

2 0.50000 0.40744 0.37357 0.35722 0.34872 0.34420 0.34177 

3 -0.38168 -0.29832 -0.26270 -0.24504 -0.23586 -0.23099 

4 0.26209 0.20479 0.17697 0.16261 0.15501 

5 -0.19443 -0.14946 -0.12729 -0.11582 

6 0.13797 0.10601 0.08972 

7 -0.10088 -0.07720 

8 0.07251 
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TABLE VI- Extrapolated values for the correlation functions Cm for the same set of 
values of A^ using the Hamer-Barber (H-B), Pade-Shanks (P-S) and Aitken-Shanks (A-S) 
transformations eqs. (9,10). 

m H-B P-S A-S 

1 -0.69829 -0.69828 -0.69827 

2 0.33890 0.33887 0.33893 

3 -0.22531 -0.22559 -0.22546 
-0.22533 -0.22536 -0.22534 

4 0.14616 0.14624 0.14620 

5 -0.10573 -0.10573 -0.10573 
-0.10352 -0.10352 -0.10352 

6 0.07279 0.07279 0.07279 
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Fig. 1- Gap and gap per spin for both models, Heisenberg (H) and Haldane-Shastry 
(HS), as a function of A^. 

Fig. 2- Correlation functions for the Heisenberg (H) and Haldane-Shastry (HS) models 
for a) \Ci\, b) C2 as a function of A^ and c) Cm for A^ = 16 vs. m. 

Fig. 3- The logarithm of the correlation function Cn/2 as a function of N for the 
Heisenberg (H) and Haldane-Shastry (HS) models. The solid line is the linear fit showing 
the exponential behavior. 

Fig. 4- a) The logarithm of the correlation functions |C(m)| andb) F{m) = ln{\C{m)\m^^'^) 
as a function of m for A^ = 16 for the Heisenberg (H) and Haldane-Shastry (HS) models. 
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